Solitons in Carbon Nanotubes 
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The symmetries of spontaneous lattice distortions in carbon nanotubes are investigated. When 
the degeneracy of the ground states remains discrete, there are solitons or domain walls connecting 
the different symmetry broken vaccua. These solitons, similarly to the case of polyacetelene, are 
fractionally charged states. In addition to the topological domain walls, there are polaron states 
with discrete energies within the energy gap. The energies and shapes of these localized mid-gap 
states should be accessible via STM spectroscopy. 
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The electronic properties of carbon nanotubes have re- 
cently become the subject of much attention B. Single 
wall nanotubes, in particular, provide a clean realization 
of quantum wires, as well as the opportunity to both engi- 
neer electronic device properties and study fundamental 
questions in low dimensional physics. Even within an in- 
dependent electron approximation, the properties of the 
single wall tubes are rather rich and useful. Depending 
on how a graphene sheet is wrapped so as to make the 
tube, the system can be either an insulator or a metal [Q. 

A great deal of work has been done forwards un- 
derstanding the role of electron-electron interactions in 
the tubes. Studies have been carried out by using the 
bosonization scheme as well as by mapping the 

problem with short range interactions into a two- leg Hub- 
bard model 0]. Within the bosonization studies, power 
law correlations were found for the order parameter of 
different electronic instabilities, such as charge density 
wave (CDW), spin density wave (SDW) and supercon- 
ductivity (SC). 

A different perspective is explored in this paper. Here 
we consider the effect of lattice deformations in the elec- 
tronic properties of the nanotubes. The effects of stress 
induced long wavelength distortions have been elegantly 
studied by Kane and Meele using a tight binding model 
f|. Also, Peierls like distortions have been previously 
investigated by studying displacements along the bond 
directions, and by assuming no spatial or quantum fluc- 
tuations of these distortions H . In this paper we show 
what exactly the symmetries of the displacement order 
parameters are, study the role of fluctuations in this order 
parameter, and discuss topological solitons and polarons 
in carbon nanotubes. 

Let us motivate the study of the symmetries of lat- 
tice distorsions by raising a question: are there frac- 
tionally charged solitons in carbon nanotubes? If there 
is a discrete number of degenerate ground states corre- 
sponding to different lattice distortions, topological ex- 
citations should exist connecting the degenerate vaccua, 
and fractionally charged excitations should be present 
in the domain walls. These ideas are familiar from an- 
other one-dimentional carbon based system: polyacety- 
lcnc |iq]. The carbon nanotubes are structurally more 



complex than polyacetelene, and this complexity will be 
reflected in the nature of the dimerization patterns that 
arise from breaking the lattice symmetries. Consider for 
example the patterns shown in Figure |l| for armchair nan- 
otubes (the tube axis is on the horizontal line). 
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FIG. 1. Dimerization patterns for armchair nanotubes. 

The ALT structures have dimerizations similar to poly- 
acetelene. However, pairs of rows along the axis have 
displacements in the opposite direction. There is a Z2 
symmetry, and there should be, at the domain walls 
between the two-fold vaccua, quantum states with frac- 
tional charge ±e/2 per spin degree of freedom for each 
of the two species of Dirac fermions present in the prob- 
lem The armchair tubes have, in the low energy 
spectrum, two species of Dirac fermions, while polyacete- 
lene has just one. The vaccum flow of charge for the 
two species of Dirac fermions goes in opposite directions, 
and it can be interpreted as flowing charge e/2 from one 
species to the other. (The precise discussion on the frac- 
tional charge will be done later in the paper). 

The Kekule bond-alternated structure contains short 
(double line) and long bonds (single line) between neigh- 
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boring carbon atoms. The patterned cells are labelled 
A,B and C according to the relative position of the sin- 
gle and double lines in the hexagons. If one visualizes the 
Kekule structure as a tiling or coloring of the hexagons 
in the nanotube with the three labels or colors A,B and 
C, there should be three degenerate vaccua correspond- 
ing to permutations of the coloring scheme. The domain 
walls between these three vaccua should be described by 
topological solitons with fractional charge ±e/3 per spin 
degree of freedom. 

The conclusions above are based on displacements only 
along the bonds. In order to fully understand the validity 
of these naive arguments, we need to look more carefully 
at more general lattice distortions and their actual sym- 
metries. We will show that the Kekule distorsion actually 
has a continuous U(l) symmetry, and fluctuations, both 
thermal and quantum, play an extremely important role. 
The ALT structure, on the other hand, is trully two-fold 
symmetric. 

Let us start from a tight-binding Hamiltonian for a 
graphite sheet: 



the cubic roots of unit z,- 



and 



zj = e lK -' Tj = e i2 3 !) . Is is also convenient to think 
of Tj as complex numbers tj = —i d Zj. Notice that 

Z)j=i z 3 = Z)j=i z j = and Sj=i z j = 3 - In tnis nota- 
tion one can write B 
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FIG. 2. Displacement vectors for the carbon atoms, (a) 
Textured displacements in the Kekule and (b) uniform dis- 
placements in the ALT structures. 



H = - E E [* + %( r )J c i W c 2 (r + Tj) + B.c. , (1) 

rGR j=l 

where r spans the triangular lattice, and the vectors Tj 
(J = 1,2,3) connect a carbon atom to its three near- 
est neighbors in the other sublattice. The distortions of 
the lattice alter the bond lengths, and thus the hopping 
matrix elements change by Stj (r) . 

In the absence of the distortions, the spectrum is given 
by E(k) = ±t \h{k)\, where ft(k) = X^ =1 e lk ^. The 
spectrum contains two Dirac points at K± = (±|^,0J. 
The dispersion /i(k) can be linearized near the Dirac 
points, i.e., k = K± + p, so the energy near these 
points is E(p) w ±i>f|p|, with a Fermi velocity Vp = |td 
(d = a/y/3 is the distance between neighboring carbon 
atoms). The nanotubes are obtained from the graphite 
sheets by wrapping around a certain direction, identify- 
ing the lattice points (0,0) and (N, M). In the N = M 
armchair tubes, the two Dirac points K± always lie on 
the p y = sub-band. 

The Kekule distortion - The size of the unit cell is 
tripled due to the dimerizations, because the hexagons 
A,B and C become distinct. This corresponds to cou- 
pling points in the original BZ which are separated by 
G = K + — K_, such as the two Dirac points. 

Consider displacements of carbon atoms that can be 
written (in terms of the undistorted lattice positions r) 
as A r = A e~ T and B r — B e* . Complex numbers 
are used to represent the displacement vectors, and A r 
and B r are in separate sublattices (see Fig. ||a). Notice 
that the textures in the two sublattices spiral in oppo- 
site directions. The Kekule distortion triples the size 
of the unit cell, hence the displacements of the three 
neighbors to any carbon atom are related by a rota- 
tion of ±27r/3. It is useful at this point to introduce 



The change in bond length dj(r), at site r and in the 



direction of Tj is 
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Using the properties of Zj , it is simple to show that the 
expression above leads to 



6dj(r) 
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(2) 



where e = ^ is the effective lattice displacement vec- 
tor that alters bonds. The other combination, namely 
j] = A ~P ' , changes bond angles without streching them, 
only costing elastic energy without any electronic gain. 
Therefore, r\ = or A = B is chosen. 
The elastic energy per hexagon is 
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Using Eq. (g) and the properties of the cubic roots of 
unit, Zj, one easily finds 5E = 3 K d 2 \e\ 2 . This energy 
cost is independent of the phase, i.e., the direction of the 
distortion of the carbon atoms. This is consistent with a 
continuous U(l) symmetry, not a discrete Z3. Terms that 
lower the symmetry appear to higher orders in the expa- 
sion of the changes in bond length (as well as bond angle). 
The non-linearities, however, are more pronounced in the 
hopping overlaps, which are exponentially sensitive to the 
changes in distance. 

Consider a change in bond hopping which is related to 
the change in bond length by an exponential: tj(r) = 
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t e aSd j( r )/ d . Expanding to second order, and using 
dj (r) as given by Eq. (|^) , one finds 



Stj(r) 



= Xj e iG r + Xj e _lti ' r + a 2 \e\ 2 , (3) 



where Xj = ^— iae — ^e 2 ^ 



The wavevector G = K + — K mixes the two species 
of Dirac fermions. Substituting Eq. (||) into Eq. (Q), one 
obtains 
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where ^ = s ip_ s A ip'_ A ), and 



iP±,s/a(p) = [ci(K± 



p)±c 2 (K ± +p)] 



are the symmetric and antisymmetric linearized fermion 
operators near the Dirac points. The order parameter is 
A/t = -3i a e + 3a 2 e 2 /2. 

The mean field gap | A| that opens is given by | A| 2 /i 2 = 
9a 2 |e| 2 + 9a 3 i (e 3 — e 3 ), and the cubic terms in e restore 
a Z3 symmetry. Notice, however, that these terms are 
smaller than the rotational symmetric leading term by a 
factor of the order |A|/i. Minimizing the elastic and elec- 
tronic energy for the filled levels, one finds (ignoring the 

non- linear effects) |A| = w^A exp ^— ^5 2 ^j-N^j , where 

«i?A is an energy cut-off scale of the order of the band- 
width t. Using typical parameters for graphite sheets 
(t re 2.4 eV, K « 19.4 eV/A 2 , a w 3.7, and d w 1.42 A), 
one finds |A| cx t e~ 2AN , so for a (5,5) tube the gap is 
of the order IK, as previously found The anisotropy 
which restores the Z3 symmetry is a factor \A\/t lower 
than the gap scale, and it only becomes apparent at tem- 
peratures of the order of 20fiK. This is a very low scale, 
and so the symmetry for the Kekule distortion is effec- 
tively U{1). 

Even at T — quantum fluctuations can restore the 
U(l) symmetry. This can be studied using a simple rotor 
model, where the arm of the rotor is the magnitude of 
the displacement of the carbon atoms from equilibrium. 
One finds that the anisotropy is irrelevant even for small 
N tubes (the estimated N c is less than 2, smaller than 
that for realistic tubes). 

Phase fluctuations of the order parameter A(x) im- 
ply a charge accumulation AQ± = ±^ Acj) |l2], where 
A<f> is the phase twist of A (2). The accumulation due 
to twisted phases of A(x) coming from the S, A chan- 
nels have opposite signs. Notice that these continuous 
phase twists, and the accompanied charge compensation 
between the symmetric and antisymmetric channels, can 
be understood in terms of a neutral boson. This is indeed 
the same situation that emerges when nearest neighbor 
electronic interactions are included, and the system is 
in the so called CDW2 phase of Krotov, Lee and Louie 



||. The charge transfered between S/A is not quantized 
because the symmetry is a continuous U(l). 

The ALT distortion - We will show that the symmetry 
for the ALT distorsion is trully a discrete Z2 symmetry. 
The change in bond length dj in the direction of Tj is now 
the same for all lattice points (see Fig. ||b), in contrast to 
the textured structure that was treated previously. One 
has 
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IZjU — IZjU , 



(5) 



where 



A-B 
2d 



Analogously to the previous case of 



lattice distortion, the elastic energy per hexagon can be 
related to u: SE = 3 (K + Kg) d 2 \u\ 2 , where in this case 
there is an extra contribution due to changes in bond 
angle, as well as bond length (Kg is defined using d to 
convert from angle to length displacements - see Ref. (ll| 
for values in graphene). Again, this energy cost is inde- 
pendent of the phase, i.e., the direction of the distortion 
of the carbon atoms. 

The correction to the Hamiltonian due to the new hop- 
ping amplitudes, however, is not independent of the di- 
rection of the displacements. It is not necessary to keep 
the changes in bond hopping beyond lowest order. Simi- 
larly to the previous case, one can show that the Hamil- 
tonian is 



H = v F J2^(P) 



h A 
-h 



*(P) (6) 



where h/ 



_ p-A 
p-A 



A = 2ia u/d . 



The spinor & = (tp]_ tl , V»+ >2 > i> 2)) wner e ^±,1/2 
are the fermions near K± in the two sub-lattices. 

For the p y — band, the distorsion opens the largest 
gap for real values of u, in which case A = — A and 
A = Wjr|j4| = 3a t \u\. There only two vaccua, corre- 
sponding to positive or negative real u (Z2), as in poly- 
acetelene. Notice that = h_^, and the spectrum has 
positive and negative energies in pairs. If u is purely 
imaginary, i.e., if the displacement is orthogonal to the 
tube axis, then there is no gap and hence no electronic 
gain from the negative energy states. There is only elastic 
cost for imaginary u, so the minimum energy path con- 
necting the two vaccua should be like in polyacetelene: a 
real u changes sign. 

Quantum numbers for ALT domain walls - The ac- 
cumulation of fractional charge in domain walls between 
the two- fold vaccua is e /2 per spin degree of freedom, and 
the ± species contribute with opposite phase shifts (one 
may allow a small imaginary part in u to see this rela- 
tive phase), hence the two quantum states have opposite 
charge. If filled or occupied, they have charge ±e/2 and 
=pe/2 respectively. One can interprete the imbalance as 
transfer of charge e/2 from one specie of Dirac fermion 
to the other. 

In the case of polyacetelene, the presence of the two 
spins (N s — 2) masks the fractionally charged states. 
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Instead, states with quantum numbers such as charge 
e and spin S = appear in the spectrum, as a conse- 
quence of having N s — 2 spin species. In the nantotubes, 
in addition to two spin states (N s — 2) there are two 
species of fermions (Nf = 2) to begin with. Because 
Nt — N s Nf — 4, the quantum numbers of the zero en- 
ergy states in the nanotubes cannot be distinguisehd from 
those of electrons. For example, one can assemble from 
the fractionally charged states an excitation with charge 
2e and spin S = on the domain wall. 

Midgap states and STM probes - In contrast to poly- 
acetelene chains, the nanotubes can be individually laid 
on a substrate, and locally probed via STM Q. One 
would then expect that the midgap states with E = cor- 
responding to domain walls (kinks and anti-kinks) could 
be probed by tunneling of electrons from an STM tip. 
The position dependent tunneling density of states would 
probe the shape of the soliton, as well as the energy of 
the state. 

In addition to the topological zero energy states con- 
necting the two ground states, there are also polaronic 
excitations. The difference between the polarons and the 
domain walls is that the polarons correspond to deple- 
tions or dimples in the order parameter without switch- 
ing between the two ground states. 

It is very simple to obtain the energy levels for the 
electronic excitations, as well as the polaron and kink 
formation energies. We start by recognizing that the 
Hamiltonian for the ALT distorsions Eq. (0) together 
with elastic energy cost 5E = 3 (K + Kg) ct|m| 2 is sim- 
ply a static version of the Gross-Neveu model for a real u 
background field |p^|-|l7||. The energies of the electronic 
states and the formation energies are 



uj n = A cos 



\2N T 



2N T 

, E n = A sin 



2N 7 



where 1 < n < Nt — 1 = 3 for the polarons, and no = Nt 
for an infinitely separated kink-anti-kink soliton pair (no- 
tice that the topological electronic state has zero energy) . 
One would hope that both the topological and the pola- 
ronic states could be probed by STM spectroscopy, with 
both energy and spatial resolution of the solitonic states. 

Interaction effects - We would like to briefly discuss 
the effects of electron-electron interaction in the above 
conclusions. The results obtained in this paper by con- 
sidering the effects of lattice displacements should be 
complementary to the studies of the electron-electron in- 
teraction effects using bosonization (Refs. The 
interaction effects also open energy gaps, of the same or- 
der of magnitude as the one discussed in this paper (see 
ref. for example). The one major difference between 
the electron-phonon interaction we discuss in this paper, 
and the electron-electron interaction discussed elsewhere, 
is that the CDW order parameter for the ALT distortion 
has a discrete Z2 symmetry. Because of the lower sym- 
metry, we expect the soliton and polaron solutions with 
mid-gap energies discussed here. 

In conclusion, we have studied the symmetries of spon- 
taneous lattice deformations in carbon nanotubes. We 



have shown that the Kekule distortion has a continu- 
ous U(l) symmetry, contrary to a naive expectation of a 
discrete Z3. Consequently, thermal and quantum fluctu- 
ations destroy any long-range order. The ALT structure, 
however, has a discrete Z2 symmetry. We discuss the im- 
plications of topological domain walls between the two- 
fold ground states, the fractionally charged states on the 
walls, and the consequences of having Nt = N s Nf — 4 
in masking the fractional states. We obtain the energies 
of the kink states, as well as polaron states, by recog- 
nizing that the Hamiltonian for the ALT distorsion is a 
version of Gross-Neveu model with a static real back- 
ground field. We argue that the mid-gap states could be 
probed by STM spectroscopy, resolving experimentally 
both the energies and the shapes of the solitonic states. 
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